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Abstract 



ZnSe scintillating bolometers are good candidates for future Double Beta Decay searches, because of the ^^Se high Q-value and 
thanks to the possibility of alpha background rejection on the basis of the scintillation signal. In this paper we report the charac- 
teristics and the anomalies observed in an extensive study of these devices. Among them, an unexpected high emission from alpha 
particles, accompanied with an unusual pattern of the light vs. heat scatter plot. The perspectives for the application of this kind of 
detectors to search for the Neutrinoless Double Beta Decay of ^^Se are presented. 
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1. Introduction 

Neutrinoless Double Beta Decay (p/3{Qv)) is, at present, a 
unique method to investigate basic neutrino properties. The 
observation of this transition would give insight on the mass 
pattern of the three neutrino eigenstates and would asses the 
Majorana character of neutrinos |[T]. The purpose of new gen- 
eration /3y6(0v) experiments [2]|3]|4l|5]El is to reach a sensitivity 
on \{my}\ of the order of 50 meV, approaching - but probably not 
completely covering - the so called inverse hierarchy region. 
The realization of an experiment with a reasonable discovery 
potential down to the smallest \{my)\ values allowed by the in- 
verted hierarchy (~ 10 meV), is today an incredible challenge 
in which detector technology plays a key role. 

In this paper we discuss the technical results obtained with 
ZnSe cryogenic particle detectors and the perspectives of their 
application to the study of /3/3{Qv) decay of the ^^Se isotope. 

2. Scintillating bolometers 

Cryogenic particle detectors (also named bolometers) have 
been proposed since many years for the study of rare events 
such as f3/3{0v) and Dark Matter 1 7 1 . The excellent detector per- 
formance and Physics results obtained by the CUORICINO ex- 
periment |8J prove how this technique is competitive in the case 
of /3/3(0v) searches. 

Unlike other solid state devices, bolometers are not based on 
ionization but on phonon detection. A bolometer is made of 
three elements: an absorber where the energy of the particles 
is deposited and converted into phonons, a phonon sensor that 
transduces the phonon signal into an electrical pulse and a ther- 
mal link to a low temperature (below 100 mK) heat sink that 
maintains the detector at the optimal temperature. 
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The working principle of these devices can be summarized as 
follow. By interacting in the absorber, a particle produces ion- 
ization and excitation, which starts triggering a series of pro- 
cesses that lead to the conversion of almost all the deposited 
energy into phonons (i.e. lattice elementary excitations). Along 
this path it is possible that a fraction of the energy is trapped 
in long-living excitated states or is lost through the emission 
of photons that escape the absorber However, a proper choice 
of the absorber (usually a single dielectric crystal, with a high 
Debye temperature and without scintillation properties) allows 
to reduce to a minimum the importance of these channels and 
therefore ensures a high efficiency in the conversion of the de- 
posited energy into phonons. The phonon signal is detected 
through a sensor that modifies its electrical properties when col- 
lided with a phonon flux. The dimensions and characteristics of 
the absorber and of the phonon sensor define the specific per- 
formance of the bolometer as far as energy resolution and time- 
response is concerned. 

In the bolometers described in this work the sensor is a semi- 
conductor thermistor glued on the crystal, while the absorber is 
a large size dielectric crystal (whose heat capacity is expected 
to dominate over the thermistor one). In this configuration the 
detector response is extremely slow: the typical rise time is of 
the order of few ms. Also CUORICINO detectors are calorime- 
ters of this kind. They are crystals of TeOa, with a mass of 
790 g. Their energy resolution at the y6y6(0v) transition is 7 keV 
FWHM, quite similar to the performance of a Ge diode. How- 
ever, when compared with the more competitive devices used 
today for /3/3{Ov) decay investigations, it is evident that the main 
disadvantages of bolometers is the lack of an active background 
rejection (e.g. particle identification as used in NEMO |9| or 
single-site vs. multi-site events distinction as in GERDA and 
Majorana |l2l|4|). The use of scintillating crystals, as absorbers, 
can overcome this limitation lfT0l[lTl[T2l[T3l . providing a pow- 
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Figure 1: Left panel: illustration of the operating principle of scintillating 
bolometers. The release of energy inside a scintillating crystal follows two 
channels: light production and thermal excitation. The heat is read out by a 
temperature sensor (NTD) glued on the primary crystal while the light is read 
by a second crystal (the light detector) where it is completely converted into 
heat. Right panel: the ZnSe Huge crystal in its Copper mounting structure. 
This crystal has few sensors glued on the top face, for redundancy. The light 
detector - here not visible - is mounted in his own Copper holder that fits in the 
top of the ZnSe one. 



erful tool to improve the experimental /3j3{Qv) sensitivities. 

The detector is in this case obtained by coupling a scintillat- 
ing crystal - operated as bolometer - to a light detector (fig- 
ure [TJ. When a particle traverses the scintillating crystal a 
(large) fraction of the deposited energy is converted into heat 
(following the chain described above), while a (small) fraction 
of it is spent to produce scintillation light (i.e. photons irradi- 
ated outside the crystal). For each interaction, two signals are 
recorded: a thermal signal produced in the scintillating bolome- 
ter and a photon signal read-out by the light detector. The ratio 
between the two (heat vs. light) depends on the particle Light 
Yield (LY). Betas and gammas have the same LY, which is typ- 
ically different from the LY of alphas or from that of neutrons. 
Consequently, the contemporary read out of the heat and light 
signals allows particle discrimination. 

In particular, if the scintillating crystal is made of a /3/3 can- 
didate, the /3/3(Qv) signal (i.e. the energy deposition produced 
by the two electrons emitted in the decay) can be distinguished 
from an alpha signal. This makes feasible the rejection of al- 
pha induced background, opening a new window in the future 
of y6y6(0v) with bolometers |fT3l . It is infact true that, in the re- 
cently concluded p/3 bolometric experiment CUORICINO, the 
major source of background is identified in alpha particles due 
to surface contaminations lfT4l [TSl [161 . This same source could 
be the main factor limiting the sensitivity of the - presently un- 
der construction - CUORE experiment 1 15 |. 

3. *^Se )8A0v) decay 

^-Se is a ySyS emitter with an isotopic abundance of 9.2% and 
a Q-value of (2995. 5±2. 7) keV ifTTl . It has always been con- 
sidered a good candidate for /3j3(Qv) studies because of its high 
transition energy and the favorable nuclear factor of merit. In- 
deed, the more recent evaluations of the Se nuclear matrix el- 
ement (NME), yield for this isotope a half-life in the range 1- 



8 lO^*- y for an \{my)\ value of 50 meV |[l8l[l9l|20l|2Tl (in the 
same condition the predicted half-life for ^^Ge is about 4 times 
larger). 

In recent years the NEMO collaboration measured its half 
life for Two Neutrino Double Beta Decay (J3/3{2v)) to be 
9.6+0.3(stat)±1.0(syst) 10''^ y and set an upper bound on the 
j6j6(0v) half-hfe of 1.0 lO^^ y at 90% C.L. |22|. 

In view of the realization of a high sensitivity /3/3{0v) experi- 
ment, the low isotopic abundance of ^'^Se is a potential problem. 
Indeed, the experimental sensitivity depends on the number of 
/3/3 emitters. With natural Se a huge mass of material would be 
needed to have a reasonable number of /3/3 emitters. However, 
this problem can be solved by isotopic enrichment, as done - 
for this same isotope - by the NEMO collaboration 191. 

On the other hand, the high Q-value is an important factor of 
merit since places the /5/3{Qv) signal in an energy region scarcely 
populated by the emissions due to natural radioactivity. This 
means that, in a calorimetric experiment, where the sum energy 
of the two electrons is recorded, a reduced number of sources 
contributes to the background and limits the experimental sen- 
sitivity |13|. The measured /3/3(2v) half-life allows to evaluate 
the irreducible background due to this decay channel, that - for 
a calorimetric experiment with high resolution - results com- 
pletely negligible (the background issue will be discussed in 
detail in section [6]l. 

4. Experimental details 

The results here discussed are a synthesis of an R&D work 
performed over more than 2 years during which we have ex- 
tensively studied the performances of ZnSe detectors in dif- 
ferent conditions. For our studies we have used three differ- 
ent ZnSe crystals whose characteristics are reported in Table [T| 
The crystals have been produced by Alkor Technologies (Saint- 
Petersbourg, Russia) and are undoped (usually ZnSe is doped 
with Te), their color variation is not fully understood, could be 
ascribed to a different stoichiometry ll23l . 

ZnSe is a semiconductor with a Debye temperature of about 
270 K [24 1 (not much different from TeOi), known as a lumi- 
nescent crystal since many years [25 1. Its X-ray induced lumi- 
nescence has been recently measured from room temperature 
down to 10 K, resulting in a dominant emission in the red re- 
gion (610 and 645 nm) and a smaller component in the infrared 
region (970 nm) [23|. The Light Yield of this crystal (mea- 
sured on small pieces cut from the same ingot from which the 
Small ZnSe crystal was obtained) was observed to increase at 
low temperature, with an emission characterized by more than 
one decay constant [26 1 . 

In all the runs discussed here, the crystals were mounted in a 
Copper frame and tightly kept in position by PTFE tips. They 
were surrounded (without being in thermal contact) by a re- 
flecting sheet (3M VM2002), but for the side which is faced 
to the light detector (LD), as shown figure [T| The tests have 
been performed operating the ZnSe scintillating bolometers at 
about 15 mK, in a low temperature refrigerator installed un- 
derground in Hall C of Laboratori Nazionali del Gran Sasso 
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Table 1: ZnSe crystal characteristics. In the last column the surface quality is 
indicated for the lateral surface and the two basis of the cylinder respectively. 
O stands for optical and R for rough. 
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quality 


Small 


yellow 


37.5 


20 


21 


ooo 


Large 


red 


120 


41 


17 


ORR 


Huge 


orange 


337 


40 


50 


OOO 



(LNGS), L'Aquila - Italy. A group of more than 10 cryogenic 
runs, with a duration ranging from 1 month to few months, have 
been realized in order to come to the results described in this 
paper In all these runs the bolometers worked properly, prov- 
ing that this crystal bears thermal cycles without a deterioration 
of its performances, whereas it exhibits an abnormal cooling 
slowness. 

Details of the set-up and the read-out are discussed in the 
next two sub-sections. 

4.1. ZnSe bolometers and the heat signal 

Each ZnSe bolometer was equipped with a Neutron Trans- 
mutation Doped Ge thermistor (NTD) [27], glued on the crys- 
tal surface and used as a thermometer to measure the heat signal 
produced by particles traversing the ZnSe. A silicon resistance, 
glued on the crystals, was used to produce a calibrated heat 
pulse in order to monitor the thermal gain of the bolometer. 
This is indeed subject to variation upon temperature drifts of 
the cryostat that can spoil the energy resolution. In most cases 
this temperature drift can be re-corrected off-line on the basis 
of the measured thermal gain variation Il28l |29l . 

The NTD thermistor transduces the heat pulse into a voltage 
signal which is amplified and fed into a 16 bit ADC. The entire 
waveform {raw pulse) is sampled, digitized and acquired. At 
the ADC input, a low pass Bessel filter is used as antialiasing to 
avoid spurious contributions in the sampled signal [30J. 

The amplitude V'"""' and the shape of the voltage pulse is then 
determined by the off-line analysis that makes use of the Opti- 
mal Filter technique ll30l[3Tl . The signal amplitudes are com- 
puted as the maximum of the optimally filtered pulse, while the 
signal shape is evaluated on the basis of four different parame- 
ters: Trise and Tdecay, TVL and TVR. 

Trise (the rise time) and Tjecay (the decay time) are evaluated 
on the raw pulse as (t9o%-tio%) and (t3()%-t9o%) respectively. For 
these detectors they are of the order of 5-10 ms and 20-30 ms. 
The rise time is mainly dominated by the time constant of the 
absorber-glue-sensor interface, while the decay time is deter- 
mined by the crystal heat capacity and by its thermal conduc- 
tance toward the heat sink. 

TVR (Test Value Right) and TVL (Test Value Left) are 
computed on the optimally filtered pulse as the least square 
differences with respect to the filtered response functionH of 



The response function of the detector, i.e. the shape of a pulse in absence 
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Figure 2: Non linearity behaviour of the ZnSe bolometer. On the X axis the 
nominal energy (E) of gamma (^^^Th and ''"K) and alpha (source A) lines. 
On the Y axis the ratio between v'"""' and E (black squares) before the en- 
ergy linearization of the bolometer response (the ratio is normalized at 1 for 
the 2615 keV peak). For the same lines the ratio between Heat and E is re- 
ported (colored circles). These data refer to the Large ZnSe crystal. A similar 
behaviour is observed for all the other detectors: the gamma pulses show a non 
linearity of the order of 3% on a 3 MeV range, the alpha lines have a similar 
non linearity but a completely different energy calibration. The experimental 
errors are of the order of 0.1-0.2% and are not visible in the plot. 

the detector: TVR on the right and TVL on the left side of 
the optimally filtered pulse maximum. These two parameters 
do not have a direct physical meaning, however they are very 
much sensitive (even in noisy conditions) to any diff'erence 
between the shape of the analyzed pulse and the response 
function. Consequently, they are used either to reject fake 
triggered signals (e.g. spikes) or to identify variations in the 
pulse shape with respect to the response function (and this will 
be our case, as in section [53] l. 

In each run, the ZnSe bolometers were calibrated with a 
^^^Th source placed between the cryostat and its external lead 
shield. Gamma lines between 238 keV and 2615 keV were 
used for amplitude (V'"""') vs. energy {E) linearization of the 
heat signal. The response of the bolometer is indeed slightly 
non linear, mainly because of the exponential dependence of 
the NTD thermistor resistance from temperature. In figure|2]we 
report the ratio V'^"" jE vs. E (where E is the full energy of the 
corresponding gamma peak) observed for one of our detectors. 
The points with £'<3 MeV correspond to gamma peaks identi- 
fied in the heat signal spectrum, the others correspond to alpha 
lines (these are due to a ^^"^Ra source facing the detector, see 
section|5]and figure|4]for more details). The non linearity mea- 
sured on gamma lines is here of about 3% on a 3 MeV energy 
range. Non linearities of the same order of magnitude are ob- 
served for all the other detectors. As it is usually done for other 
bolometers |fT6ll32l . the non-linearity is corrected by introduc- 
ing the calibrated amplitude of the heat pulse {Heat) - obtained 



of noise, is computed with a proper average of a large number of raw pulses. It 
is also used, together with the measured noise power spectrum, to construct the 
transfer function of the Optimal Filter. 
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Figure 3: Th gamma calibration of the Huge ZnSe crystal. The filled his- 
togram is the projection - on the Heat axis - of the events visible in the scatter 
plot of figure|5](right panel). The black line histogram is obtained as projection 
on the Heatg axis (wher e 6 ha s the value that optimize the RMS of the 2615 keV 
^"^Tl peak), see section [sT] for more details. In the inset the highlight of the 
2615 keV line. The gain in resolution on all the lines is evident. 



from V - using the following relationship: 

log(Heat) = A + B log(Vh,at) + C logiVhemf 



(1) 



where the three coefficients A, B and C are determined fitting 
the (V*'''", Ei) points of the calibration gamma lines. It is worth 
to note that in this way we attribute to the heat signal the nomi- 
nal energy E of the detected gamma line, i.e. the cahbration 
returns Heat-E for gamma particles. Here E is the total energy 
lost within the crystal, while only a fraction of it is converted 
into heat. The other fraction is spent in scintillation light or in 
possible blind channels such as long living metastable states. 
The result obtained after the calibration is shown in figure |2j 
the colored circles are the ratio HeatjE for each gamma and 
alpha line. The calibration returns, as expected, Heat/E~l for 
all the gamma lines. While for alpha lines Heat/E >1 (all the 
points with E > 3 MeV in figure[2|. This means that alpha lines 
appears in the Heat spectrum with a displacement of about 1 
MeV, i.e. for the same total energy deposition in the crystal, 
heat pulses due to alphas are higher (by a factor ~1.15) than 
heat pulses due to gammas. Provided that we observe for alpha 
particles a higher light signal, this behaviour is quite strange. 
We will come back to this feature in section[ 



only used it at the beginning of our studies on ZnSe, we did not 
have the occasion for a long calibration measurement. 

4.2. Light detectors and the photon signal 

The scintillation hght emitted by the ZnSe crystals is de- 
tected with small bolometers ||33l . These are made by high 
purity Ge wafers with a large area (to optimize photon collec- 
tion) but a small mass (to reduce the heat capacity and therefore 
to have a sizable temperature increase upon photon absorption). 
They are 36 mm in diameter and 1 mm in thickness. The sur- 
face of the Ge wafers facing the scintillating crystal are dark- 
ened through the deposition of a 600 A layer of Si02 in order 
to increase the light absorption. We have used also a very large 
area Ge wafer (66 mm in diameter) able to cover contemporary 
the Small and Large ZnSe crystals in order to provide a direct 
comparison of their Light Yields. 

The Ge wafers are provided, as the ZnSe bolometer, with 
an NTD thermistor for the temperature read-out and with a Si 
resistor for gain monitoring. Electronics and DAQ are identical 
to that of ZnSe scintillators. Rise and decay times of the signal 
are of about 2 and 15 ms respectively. 

The Ge wafer are too small to be able to detect the ener- 
getic gamma lines of the ^^'^Th source used for ZnSe calibra- 
tion. Therefore a weak ^^Fe source was faced to the Ge wafers, 
on the side opposite to the crystal scintillator, for calibration 
purposes. The 5.9 keV Kq, X ray line of ^^Fe was used to mea- 
sure the energy resolution of the detector and its amplitude to 
energy conversion. Obviously, with only one calibration peak, 
we are not able to correct the Light amplitude for the non linear- 
ity of the LD response. However, the expected non linearity is 
similar to the one recorded for the scintillating bolometer (since 
the NTD sensor is similar) while the energy range is smaller by 
more than one order of magnitude. The effect should be there- 
fore negligible. 

The ^^Fe calibration provides a method to measure and com- 
pare the LY's of different crystals. This cannot be used for an 
absolute evaluation of the amount of scintillation light emitted 
by the crystal since the Light amplitude is not corrected for the 
photon collection efficiency that is certainly below 1. Finally 
the energy resolution of the LD is measured as the FWHM 
of the 5.9 keV line. Depending on the detector used and on 
the working conditions, the resolutions range from 200 eV to 
600 eV FWHM. These resolutions are negligible if compared 
with the width of gamma and alpha lines in the Light spectrum. 



The energy resolution on the heat signal was evaluated mea- 
suring the FWHM of the 2615 keV line in the Heat spectrum. 
The peak FWHM was evaluated by a fit with an asymmet- 
ric gaussian. The best resolutions measured (average of left 
and right FWHM) for the Large and Huge crystals were of 
21 ±2 keV and 28 ±1 keV respectively (see filled histogram in 
figure |3]l. 

In the case of the ZnSe Small crystal we have not been able 
to evaluate the energy resolution because no gamma peak is vis- 
ible in the calibration or background spectra. This is explained 
by the small dimensions of the crystal and its low Z. Since we 



5. Results 

We discuss in this section the results obtained in several runs 
performed in the last few years, with the three above mentioned 
crystals. Our measurements aimed mainly at the study of the 
alpha vs. beta/gamma discrimination properties of ZnSe, at the 
evaluation of the LY and of the Quenching Factor (QF), and at 
the study of the variation of these properties upon variation of 
the crystal color and surface quality. In order to be able to com- 
pare the detector response to alpha and beta/gamma particles, in 
most of the runs, the crystals were faced to low intensity alpha 
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Figure 4: Spectmm of source A as recorded by a Si surface barrier detector 
(having 35 keV FWHM energy resolution at 5 MeV). The energy [keV] of the 
main lines is indicated. 



sources (while gammas were generated with the external ^^^Th 
source used for calibration). Two different radio-isotopes have 
been used so far. 

Source A was obtained by implantation of ^^"*Ra in an Al re- 
flecting stripe. The implantation was accomplished by exposing 
the stripe to a thin ^^**Ra source. The ^^'^Th isotopes contained 
in the source alpha decay to ^^^Ra. The recoil energy of ^^"^Ra 
is sometimes enough to eject the nucleus from the source and 
implant it in the Al stripe. The extremely short range of nuclear 
recoils in solid materials allows to obtain a very shallow depth 
source (this can be appreciated in figure|4]from the symmetry of 
the peaks). Besides several alpha lines of different relative in- 
tensity, this source also produce beta particles with a maximum 
energy of 5 MeV due to the beta decay of ^'"^Tl. 

Source B was obtained contaminating an Al stripe with 
a ^^^U liquid solution and later covering the stripe with an 
alluminated Mylar film (6 fxm thick). Thus the source has a 
continuous spectrum extending from about 3 MeV to 0. 



The usual way to visualize alpha vs. beta discrimination is to 
draw the Light vs. Heat scatter plot. Here each event is identi- 
fied by a point whose abscissa is equal to the Heat signal am- 
plitude (recorded by the scintillating bolometer), and ordinate 
equal to the Light signal amplitude (simultaneously recorded 
by the LD). In the scatter plot beta/gamma, alphas and neutrons 
give rise to separate bands, in virtue of their characteristic Light 
to Heat ratio. As already anticipated, this feature is the result 
of the different LY's characterizing these particles. 

The Light vs. Heat scatter plot of the ZnSe Huge crystal, 
exposed to the calibration gamma source (^^^Th) and to the al- 
pha source B, is shown in figure |5] This scatter plot (as well 
as all the others we obtained with ZnSe crystals) is completely 
different from any other we ever measured with different scin- 
tillating crystals (we tested several scintillating bolometers like 
CdW04, ZnMo04 and CaMo04). Three anomalies are indeed 
here evident: 



1. monocromatich gamma peaks lie on positive slope lines 
(right panel of figure |5]), whereas a negative slope is ex- 
pected; 

2. the alpha band lies above the beta/gamma band, suggesting 
an alpha QF larger than one (left panel of figure|5]); 

3. alphas are always accompanied with a tail that draws neg- 
ative slope lines (figure |8]l. 

The anomalies resulted reproducible for the three crystals, 
despite the different crystal colour and size and despite the dif- 
ferent experimental conditions. In the following three sections 
we document these unexpected features and try to give a possi- 
ble interpretation of them. 

5.7. Correlation between Heat and Light and energy resolution 
on gammas 

In a standard scintillator, monochromatic events (i.e. those 
corresponding to the complete absorption of a monochromatic 
particle in the scintillating crystal) appear as sections of lines 
with negative slope, showing an anticorrelation between Heat 
and Light Il32l[34l . The observed spread along a negative slope 
line can be accounted for assuming the existence of fluctua- 
tions in the Heat and Light signal amplitudes. Since the energy 
has to be conserved, the amplitude of the fluctuations (SLight 
and 5Heat) must compensate each other: 6Light--6Heat. In 
other words the fluctuations are anti-correlated between each 
other and produce the observed pattern. Finally, they are vis- 
ible because their amplitude is much larger than the intrinsic 
resolution of the detectors. This is what happens for example 
for CdW04 crystals. In this case a rotation of the axis can be 
used to improve the energy resolution of the detector |32, 34 1. 

In ZnSe crystals the slope of gamma lines is clearly positive 
(see right panel of figure [5]): for a monochromatic event (as it 
is the photopeak of a gamma line) the minimum heat release 
corresponds to the minimum of light emission, and vice-versa. 
It seems that fluctuations in the energy conversion into heat and 
light are correlated, and that there is a lacking energy which is 
feeding some invisible channel. 

The data presently in our hands are not enough to proceed 
further in the modeling of this unusual behaviour of ZnSe scin- 
tillator, we simply observe that - also in this configuration - a 
rotation of the axis 1 32 1 can be used to improve the energy res- 
olution of the detector. 

Figure |3]shows (filled histogram) the spectrum obtained pro- 
jecting the gamma events of figure [5] on the Heat axis. The 
black line histogram is obtained by projecting events on a ro- 
tated axis (Heatg) defined as: 



Heatg — Heat cosiO) — Light siniO) 



(2) 



where the rotation angle is chosen in order to reach the mini- 
mum FWHM for the 2615 keV. The improvement obtained on 
the resolution is of about a factor 3 on the 2615 keV line of 
^o^Tl (the FWHM is reduced from 28±1 keV to 9.5+0.4 keV). 
On the other peaks the improvement is smaller: at 911 keV from 
7.8+0.6 keV to 5.9+0.3 keV, at 583 keV from 7.5+0.5 keV to 
5.4+0.3 keV. Provided that it is proved that also /3/3(Qv) events 
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Figure 5: Scatter plot Light vs. Heat recorded with the Huge ZnSe crystal. The detector was here exposed to the B alpha source and to the calibration gamma source 
(^^^Th located outside the refrigerator). Three populations of events have been marked: black spots for beta/gamma's, colored spots for alphas, circles for events 
with an ionization component in the LD. In the right panel a zoom of the beta/gamma band. The gamma lines due to ^'^Th are clearly visible. The slope of all the 
lines is positive. 



Crystal ZnSe Small ZnSe Large ZnSe Huge 

LY [keV/MeV] 1.3 7.5 4.6 

QF„ 4.4 4.2 3 



Table 2: LY^ measured on the 2615 keV line of ^"^Tl (^^^Th) and QF„ mea- 
sured for the 5.7 MeV alpha line of ^-'*Ra. For each crystal we report values 
measured in similar experimental conditions: with the LD mounted on the top 
of the crystal and the alpha source facing the opposite side. The experimental 
errors on LYy and QFn, are of the order of 5% and 10% respectively. In the case 
of the ZnSe Small and Large crystals we were also able to compare directly 
their LYy in a dedicated run, where their scintillation photons were detected by 
the same Ge wafer (6 mm in diameter). The results where consistent with what 
reported in this table. 

would give rise to a monochromatic line with the same slope, 
this technique could be used to improve the resolution of the 
detector. 



particles of the same energies. It is however difficult to have al- 
pha and beta/gamma sources with monochromatic lines of the 
same energies therefore we have evaluated the QFq, as: 

^ Ught^ Energy (4) 
Energya Lighty 

where we always compare the 5.7 MeV alpha line of the im- 
planted source with the 2615 keV gamma line of the ^^^Th 
source. 

The QF obtained for the Small, Large and Huge crystals, for 
an identical position of the alpha source (i.e. facing the crystal 
basis opposite to the LD) are reported in table[2] QF„ values are 
definitely larger than 1, and this unusual feature is accompanied 
by other peculiar behaviours that we extensively investigated in 
several dedicated runs. The experimental results and the con- 
clusions that can be drawn are discussed in the next subsections. 



5.2. Light Yield and anomalous light emission from the alphas 

The Light Yield of a crystal is conventionally defined as the 
fraction of the total energy deposited in the scintillator which is 
converted into light. As discussed previously, we do not have 
an absolute calibration of the Light signal, therefore we can 
provide only a relative evaluation of the Light Yield defined as: 



Lighty 
Energjy 



(3) 



A similar definition holds for alphas. In our case, we chose 
to evaluate the LY^ on the 2615 keV fine of ^"'^Tl (i.e. 
Energyy=2.615 keV). The results for our 3 crystals are reported 
in table 121 

The QF„ is usually defined as the ratio between the ampli- 
tudes of the light pulses corresponding to alpha and beta/gamma 



5.2. L The anomalous alpha QF 

To our knowledge no other scintillator with a QFq. definitely 
higher than 1 have ever been reported in literature, and this 
feature cannot be easily accommodated within the theoretical 
framework used to describe the scintillation properties of ma- 
terials. Indeed heavy particles usually have QF lower than one 
and this is explained on the basis of a saturation effect due to 
the high ionization density that characterizes the interaction of 
heavy particle in matter 1351. 

We have considered the possibility that a QFq. larger than 1 
could be due to the experimental configuration, in some way 
privileging the collection of photons emitted by the point-like 
external alpha source with respect to the light produced in crys- 
tal bulk by gamma rays. We investigated therefore the light 
collection, the self-absorption of the scintillation photons in- 
side the crystal and - finally - the transparency of the LD's to 
the scintillation photons. In more detail: 
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• a self-absorption of the light inside the scintillating crys- 
tal could penalize gamma interactions (that take place in 
crystal bulk) with respect to alpha ones (alphas come from 
external sources and interact in a thin surface layer). To in- 
vestigate this effect we placed the alpha source on the side 
opposite to the LD in such a way that only the light travers- 
ing the crystal could be detected. All the QF's reported in 
table |2]have been measured in this configuration. 

• to check possible problems arising from light collection 
we compared a measurement done with a ZnSe Huge crys- 
tal and a CdW04 crystal of the same size. Both crystals 
were mounted in a completely identical configuration with 
the same LD. The result was a normal behaviour in the 
CdW04 crystal and an anomalous one in the ZnSe. The 
LYy measured for the CdW04 crystal was 17 keV/MeV 
||32]| . Consequently, ZnSe has a light emission on gam- 
mas that is four times lower than the one recorded with 
CdW04, while the light emission on alphas is nearly two 
times larger. 

• we investigated the possibility that the low LYy we mea- 
sure could be due to a dominant beta/gamma photon emis- 
sion on a wavelength to which our LD is partially transpar- 
ent. The Ge wafers are completely opaque to the red light 
that, as measured upon X-rays excitation |^3l, appears to 
dominate the ZnSe luminescence spectra. Nonetheless, at 
the very low temperature where we work, something could 
change. To verify this hypotesis we made a test using - as 
LD - a sapphire bolometer covered with a lead substrate 
(completely opaque from IR to UV). No relevant modifi- 
cation in the measured LY's and QF's was observed. 

All our measurements, and consequently the above reported 
considerations, have been based on external alpha sources since 
we have never observed - in our crystals - an internal alpha line. 
We investigated the possibility that bulk alphas do not scintil- 
late, in this case we should see heat signals in anticoincidence 
with light ones. These dark counts are present in our spectra, 
but with the same percentage in the alpha and gamma region. 
They can in fact be ascribed to a light-heat coincidence lost due 
to trigger inefficiency on one of the two channels. 

As a final remark on this topic, it is worth to recall that the 
value of the QFq, may depend on the parameter that is used to 
measure the energy content of the recorded pulse. In our case 
this is the pulse amplitude. In a very naive model of the bolome- 
ter, that assumes an instantaneous energy deposition in a mono- 
lithic device ll36ll . the amplitude is proportional to the energy 
released in the crystal. However, we have clear indications that 
these approximations do not hold. As it will be discussed in 
section |5.3| we observe a clear difference between the shape 
of alpha and beta/gamma light pulses. Figure [6] shows aver- 
aged alpha and beta/gamma light pulses having the same Heat 
amplitude (and consequently produced by particles releasing in 
the scintillator almost the same energy)Fl Once normalized to 



-The time window is here 4 sec. wide. The average was needed to reduce 




Figure 6: Averaged light pulses from beta/gamma (black line) and alpha (col- 
ored line), after pulse height normalization. In shadow the response function of 
the LD, as obtained averaging pulses due to direct ionization in the Ge wafer. 
The real shape of the light signals shall be obtained by deconvolving the signal 
with the response function of the detector. The time window is 4 sec. wide. 

the same Light amplitude the two signals show a quite differ- 
ent area. If the area, and not the amplitude of the pulse, would 
be used as an estimator of its energy content the QFq, would 
result lower than reported in table |2] (by about a factor 1.5), al- 
though still much higher than I. It could therefore be argued 
that, upon collection of all the light emitted by beta/gammas, 
the QF„ could approach a more standard value. 

For the sake of completeness, we mention here that a 
strange behaviour of the QF„ in ZnSe:Te (Tellurium doped 
Zinc Selenide) was already reported in literature. In particular, 
E. V. Sysoeva et al. JST) quote a QF„=1.05. While W. Klamra 
et al. fSS] quote a QFq.=0.74-0.69, depending on the shaping 
factor of the amplifier. It is however very difficult to compare 
our results with those reported in literature. There are indeed 
two key aspects that have to be pointed out: 

• all the measurement performed before this work were car- 
ried on extremely thin crystals (1-2 mm). This because the 
transmission band of ZnSe - at room temperature - is very 
close to the emission band (while at low temperatures the 
crystal becomes almost transparent at the 610 nm wave- 
length of the dominant emitted radiation). 

• our measurements are done on undoped ZnSe crystals at 
extremely low temperatures, while most of the literature 
concentrates on doped samples measured at much higher 
temperatures. In our case the choice is motivated by the 
fact that the doping introduces a strong modification of the 
crystal lattice, resulting in a larger heat capacitance, thus 
strongly decreasing the performances of a bolometer On 
the other hand the use of an undoped crystal is advanta- 
geous because - as well established - decreasing the tem- 



the noise contribution and clearly appreciate the long tail of the pulses. It has 
to be remarked that - in the amplitude interval where we averaged pulses - no 
dependence of pulse shape on amplitude is observed (i.e. the average does not 
change the shape of the signal). 
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Figure 7: Scatter plot of Light vs. Heat recorded with the Huge ZnSe crystal in 
two different runs. In run I (colored points) the B alpha source was mounted on 
the bottom side of the crystal, opposite to the LD. The scatter plot obtained in 
this run is reported in figure[5] In run II (black points) source B was moved on 
the side of the crystal and in its place was mounted source A. The dotted lines 
identify the opening of the alpha band on a wide angle, due to the variation of 
the QF„ upon source position. 



perature the LY of intrinsic scintillators increases, while 
the opposite happens for doped scintillators. 

5.2.2. QF dependence on source position 

Figure |7] compares the results obtained with the same detec- 
tor changing only the position of the alpha sources. In run I 
the B alpha source was mounted on the bottom side of the crys- 
tal, opposite to the LD (this is the same scatter plot reported in 
figure[5]l. In run II source B was moved on the side of the crys- 
tal, while in its previous position was mounted source A. By 
comparing the QFa^ obtained for the two positions, it is clearly 
evident a variation of the order of 20% . This indicates that - 
in the case of a source distributed uniformly around the crystal 
(and not point-like sources like those we are using)- we should 
observe an alpha band spanning a much larger angle. 

In figure |5] and figure |7] we observe the presence of back- 
ground alpha events forming a positive slope line (similar to 
the ones observed for gamma particles, see section 5.1 1 at 



Heat~6.2 MeV and Light~^0 keV. This region corresponds to 
a nominal alpha energy of about 5.3 MeV, which is the en ergy 



4.1 



of the alpha decay of ^"'Po. Indeed, as discussed in section 
the correct energy of an alpha particle is lower than the recorded 
Heat amplitude by about a factor 1.15. This line can be as- 
cribed to a ^"'Pb contamination (a quite common contaminant 
whose strongest signature is the alpha emitted by ^"^Po ) be- 
cause its intensity does not change in time (^"'Po has a half- 
life of 138 days and we have two measurements with the same 
crystal at a distance of ~4 months, without sizable variation of 
intensity). The line has an evident tail (with negative slope) to- 
ward low Light - high Heat signals, an indication (as wiU be 
discussed in section [5.2. 3| l that it is due to a surface contami- 
nation. A reasonable assumption is that the source is a ^"^Pb 
surface contamination of the crystal or of the mounting struc- 
ture surrounding the crystal. In this case, the extension of the 



line on the Light axis can be interpreted as due to the variation 
of the QF„ as a consequence of a position dependence. In fig- 
ure |7] we have indicated approximately the width that we could 
expect for the alpha band assuming that the ^'°Po line is due to 
an uniform contamination surrounding the crystal. 

The dependence of QF„ on the source position could in prin- 
ciple be due to a difference in photon collection efficiency, how- 
ever in a test run done with two LD (where therefore the col- 
lection efficiency was definitely improved) the extension of the 
^'°Pb line on the Light axis was not observed to change. More- 
over, upon changes in the alpha source position, we observe a 
variation not only in QF„ but also in pulse shape (the decay 
time of the light signals). All this gives credit to a second hy- 
pothesis: that the variation of QF„ upon source position is due 
to position dependence of LY„ (a similar effect is discussed in 
]) or to the self-absorption of the emitted light. 



5.2.3. Low Light tails of the alphas in the scatter plot 

As anticipated in the previous section, the separation of the 
alpha band from the beta/gamma one is not complete. This is 
evident looking at monochromatic alpha lines, like those pro- 
duced by source A. For each alpha line, we observe a spot plus 
a long tail of events - with lower Light and larger Heat - that 
draws a negative slope line on the scatter plot (see figure |8]l. 
Obviously the same effect produces - in the case of source B - 
a continuum without evident structures below the alpha band. 
This feature is potentially dangerous since the purpose of our 
search is the development of a detector able to identify and re- 
ject exactly surface alphas. Indeed, this is the most pernicious 
background source limiting the sensitivity of the CUORICINO 
experiment. We have observed that this problem appears to be 
connected to the quality of crystal surface (or to the presence - 
on the surfaces - of residuals of the abrasive, non-scintillating, 
powder used for polishing). Indeed the fraction of events pop- 
ulating the tail toward low light emission is different (for the 
same crystal) if the alpha source hits an optical face of the de- 
tector or a rough face. Figure|8]illustrates the test we have done 
and its results. The problem persists, although in reduced size, 
also in optical crystals as is the case of the ZnSe Huge. 

5.3. Pulse shape analysis and alpha particle discrimination 



As discussed in sections 4.1 and 4.2 the time response of 
the bolometers (both the LD and the scintillator crystal), is ex- 
tremely slow. The light pulse has rise time of the order of few 
ms that - compared with the typical response of a photomulti- 
plier or a photodiode - is extremely slow. However, at low tem- 
perature, also the light emission from scintillators is character- 
ized by long time constants, with different values for alphas or 
beta/gamma's [|40l . This happens also in ZnSe (as evident from 
figure |6]l and a discrimination based on pulse shape analysis is 
achievable. 

Figure |9] shows the decay constant (left panel), of events 
recorded with the LD, versus the energy of the corresponding 
heat pulse (i.e. Heat). This plot refers to the measurement done 
with the ZnSe Huge crystal, exposed to a ^^^Th gamma source 
and the B alpha source, whose scatter plot is shown in figure |5j 
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Figure 8: Scatter plot of Light vs. Heat recorded with the Large ZnSe crystal. 
The detector was here exposed to the A alpha source and to the calibration 
gamma source (~^^Th located outside the refrigerator). Each alpha peak comes 
with a - low Light high Heat - tail that draws a negative slope line in the scatter 
plot. In the inset: histogram of the events belonging to the 5.7 MeV alpha line 
(colored spots) projected along their negative slope line. The two histograms 
have been obtained for the same source (source A) and the same detector. The 
only difference being that in the case of the black histogram the non optical 
basis of the ZnSe cylinder was in front of the source, while in the case of the 
colored histogram it was the optical basis that faced the source (the scatter plot 
refers to this latter configuration). 



The three populations identified in figure [5] appear to have a 
clearly diff'erent light pulse shape: 

• events of direct ionization on the LD (circles) are very fast, 
with a decay time of about 5 ms. These events are interac- 
tion that involved both the LD and the scintillator There- 
fore they are in coincidence with a heat pulse in the scin- 
tillator, while the energy recorded by the LD is in part due 
to direct ionization and in part to scintillation light. Also 
^^Fe events have this same decay time (they are not visible 
in figure |9]because they are not in coincidence with a heat 
pulse). 

• alpha events (colored spots) have an average decay time of 
less than 10 ms and have a distribution quite different from 
that of beta/gamma events (black spots). Events belonging 
to the ^'"Po line appear to have a wider distribution of the 
decay time than events produced by source B alphas. This 
looks as the same effect that, on the scatter plot of figure]?] 
makes ^"'Po events extend on a wide range of Light am- 
plitudes. 

• beta/gamma events (black spots) have an average decay 
time of less than 13 ms with a wide distribution that is 
partially superimposed to the alpha one. 

A similar behaviour is observed on the rise time. 

Also the shape of heat pulses appears to be different for alpha 
and beta/gamma events. In the right panel of figure|9]it is shown 
the rise time of heat pulses. Beta/gammas (black spots) have a 
faster rise time than alphas (colored spots), the same is true for 
the decay time. 



The reason for having a faster light response and a slower 
heat one in the case of alphas is not clear and will require dedi- 
cated measurements and the construction of a model for signal 
formation. This will probably help also in the understanding of 
the other anomalies of this scintillator, discussed so far. 

The observed difference in pulse shapes, open in any case 
an important possibility for the use of ZnSe bolometers for 
PfiiQv) search: the rejection of alpha particle background on the 
basis of the pulse shape. Indeed, while the use of the Heat to 
Light ratio seems to provide only a partial discrimination (due 



to the low Light events discussed in 5.2.3 1, the use of the pulse 
shape of light pulses appears much more efficient. 

In a dedicated run - with optimized measuring conditions - 
we were able to reach a quite high alpha rejection efficiency. 
The ZnSe Large crystal (the one with the best LY) was mounted 
with - on its bottom - source B and a Sm source ('"^^Sm alpha 
decays with the emission of a 2.2 MeV alpha particle). The 
LD was operated at a slightly higher temperature than usual, in 
order to have a faster response. The signal sampling was in- 
creased to 16384 points on a 4 sec window, to improve the res- 
olution on the signal shape. The result is illustrated in figurefTOl 
The scatter plot (left panel) of Light vs. Heat shows clearly the 
alpha and beta bands. The alpha band has a quite important tail 
towards low Light signals and - based on the Heat/Light ratio - 
the alpha discrimination looks quite poor The LD clearly dis- 
tinguishes alphas from betas on the basis of their pulse shape 
(central panel), here measured with the TVR parameter. The 
projection of the TVR distribution (right panel) shows how ef- 
ficient can be this kind of discrimination; for a better visuali- 
sation of this efficiency all the signals with TVR lower than 2 
appear in the scatter plot (left panel) in color The alpha rejec- 
tion efficiency was evaluated in the region 2400-2600 keV (the 
nearest region to the **^Se /3/3 Q-value and - in our measurement 
- populated enought by both beta/gamma's and alpha's) to be 
about 15 sigma. 

Finally, the possibility of an alpha/beta discrimination based 
on the scintillating bolometer only (as that presented in the right 
panel of figure [9]) is extraordinarily interesting since implies an 
oversimplification for the detector (the LD is no more neces- 
sary). However, in our measurements, it is still not enougth 
sensitive and reproducible. We recall here that ZnSe is not the 
only scintillating crystal were we have observed this character- 
istics; a very promising crystal is for example CaMo04 BTIl . 

6. jS/S(Ov) perspectives with ZnSe 

In this section we briefly discuss the perspectives in the ap- 
plication of ZnSe scintillating bolometers for /3/3{Ov) study. To 
do this we assume a FWHM energy resolution of 10 keV on 
the Heat signal at the Se Q-value (-2995 keV) and a ROI 
two FWHM wide (i.e. from 2985 to 3005 keV). For the al- 
pha rejection we assume to use the pulse shape analysis on the 
light pulse, with the same efiiciency reported in section [53] (15 
sigma). 

In order to translate the results discussed here in the actual 
sensitivity of a j3/3 experiment based on this technique, we can 
refer to the only existing cryogenic setup, the one installed at 
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Figure 9: Pulse shape of the events represented in the scatter plot of figure|5] Left panel: decay time (Tjgca\) of light pulses vs. energy (Heat) of heat pulses. It 
is here evident that the three populations identified in figure|5]are characterized by different distributions of Tjgcav. Beta/gamma particle events (black spots) have 
an average rjecar of 13 ms, with a wide distribution that partially overlaps with the alpha event distribution. Note that the difference in the width of the alpha and 
beta/gamma rj^cav distributions is simply due to the dilferent weight of noise on their light pulses. Alphas have a much higher Light amplitude that beta/gamma's, 
therefore the effect of noise is smaller. This is not evident from the plot since we show Tj^cav of light pulse vs. the heat pulse amplitude, (Heat). Right panel: rise 
time (t,.;,,,) of heat pulses vs. energy (Heat) for these same events. Alpha and gammas shows different distributions, that however overlap. 




Figure 10: ZnSe Large exposed to source B and to a Sm source, both mounted at crystal bottom. Left panel: scatter plot Light vs. Heat. Black spots are alpha 
particles and colored spots beta/gamma's. The distinction is based on a cut defined on the TVR parameter on light pulses. Central panel: TVR distribution for 
light pulses vs. the amplitude of the corresponding heat pulse (Heat). The cut that is used to separate alphas from beta/gamma's is TVR=2. Right panel: TVR 
distribution for Light pulses between 2400 and 2600 keV. 



10 



LNGS (Hall A) in the 80's by the Milano group and later used 
also for the CUORICINO experiment. The background model 
for this setup is in fact well known |42 1 and the extrapolation to 
the technique described here can be obtained by simply substi- 
tuting TeOa crystals with ZnSe ones. 

For the background rate induced by environmental neutrons 
and muons we use the results reported in Ii42il and 143 J . In- 
deed the change in the molecular compound does not affect 
in a relevant way the background rates. External neutrons 
yield a contribution integrated in the 3-4 MeV region lower 
than 10"^ c/keV/kg/y ll42ll . while external muons - operating 
the array in anticoincidence - yield a contribution of about 
8 IQ-^ c/keV/kg/y in the 2-4 MeV region [43\. 

For the background induced by environmental gammas, the 
results reported in [42 1 do not consider the energy region above 
the 2615 keV line of ^"'^Tl where the counting rate is drastically 
reduced (as evident in figure 1 of |42|). In our ROI the dom- 
inant contribution comes from ^''*Bi and ^"^Tl. Indeed ^'"^Bi 
has few low intensity gamma rays in or above our ROI (the to- 
tal branching ratio of these gamma lines is 0.08%) while ^^^Tl 
has no gamma with energy larger than 2615 keV, but can con- 
tribute to the counting rate in our ROI through the contempo- 
rary emission of two photons with a sum energy of 3.2 MeV 
(the 2615 keV and 583 keV gammas). The background model 
developed for the Hall A set-up at LNGS, and described in 1 42 J . 
identifies the major source of ^'"^Bi and ^osf j activity as due to 
contaminations localized in between the external and the in- 
ternal lead shields. This allows us to evaluate ^'"^Bi and ^''^^Tl 
contributions in the ROI. These are: ~0.006 c/keV/kg/y from 
2i4Bi (dominated by the 3000 keV Une) and 3 lO^^ c/keV/kg/y 
from 2"8ti . 

Concerning the background contribution of the detector it- 
self, the 15 sigma efficiency in alpha rejection reduce - for rea- 
sonable contamination levels of the surfaces and of the crys- 
tal bulk - to absolutely negligible values the alpha rate in the 
ROI. Finally, the use of delayed coincidences allows to get 
rid of the most worrisome background arising from ^^'^U and 
■^^■^Th crystal bulk contaminations. The measured upper lim- 
its on the ZnSe (the Huge crystal) internal contaminations are 
of 0.4 ;uBq/kg (at 90% C.L.) for both the two chains; these 
contaminations results in a background counting rate far below 
10-"* c/keV/kg/y. 

In conclusion the ultimate sensitivity of a ZnSe based 
j6/S(0v) experiment performed in the Hall A set-up at LNGS 
would be limited mainly by the ^'"^Bi contamination that leads 
to a background rate in the ROI of ~0.006 c/keV/kg/y . 

7. Conclusion 

ZnSe scintillating crystals present a number of anomalies or 
unexpected behaviours that certainly need further investigation 
in the future. The QFq, larger than one could find an explana- 
tion in the existence of a light emission from gammas on ex- 
tremely long time scales (something like an after-glow). How- 
ever - even if the usual ratio between LY^ and LY^ - could be 
restored with proper collection of all the emitted light, the larger 
amplitude of alpha particle heat signals with respect to gamma 



heat signals (see figure|2]and section 4.1 1 remains unexplained. 
The same holds for the positive slope of gamma monochromatic 
lines in the scatter plot (and could be also for internal alphas). 
The understanding of these behaviours is obviously mandatory 
when projecting a /3/3 experiment since they can influence the 
response of the detector to a /3/3(Qv) decay and therefore its ef- 
ficiency and resolution. Nonetheless some conclusion can be 
drawn concerning the possibility of employing this crystal for a 
high sensitivity fifiiQv) experiment. The energy resolution - pro- 
vided that the technique of axis rotation is proven to be reliable 
also for a y6yS(0v) event - is comparable with that of CUORI- 
CINO detectors and not very far from Ge diodes. As a prove 
of the perspectives for the application of this detector to future 
/3/3(Qv) searches, we have evaluated the background achievable 
with a medium size experiment in a realistic set-up. We con- 
sidered the existing cryogenic set-up of Hall A in LNGS that 
was specially built - in the eighties - for low background exper- 
iments. The result is a counting rate of about 6 lO""* c/keV/kg/y 
in the ROI, dominated by the radioactive contamination of the 
set-up itself. However according to the measured bulk contam- 
ination of the ZnSe Huge crystal, the intrinsic background of 
the detector is lower than lO"** c/keV/kg/y, proving that - in an 
extremely low background set-up - very high sensitivies for the 
PPiOv) process are achievable. 
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